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r ~T l Abstract 

Based on a neutrino mass matrix model in which a 2-3 symmetry is only broken by 

Oh| a phase parameter, it is investigated how the lepton mixing matrix can deviate from the 

D 

fH ' so-called tribimaximal mixing under a condition that CP is maximally violated. 

^ \ 1 Introduction 

<3\ [ Recent data in neutrino oscillation experiments [TJ [2] have showed that the lepton mixing 

, matrix U is fairly in favor of the so-called tribimaximal mixing [3] Utb- If U is exactly U = Utb, 

then CP violation in the lepton sector is absolutely forbidden. Our interest is in how the CP 

violation in the lepton sector can be large and how U can deviate from the tribimaximal mixing 

G\ '. Utb- 
O ■ 



X 



Meanwhile, we know that a neutrino mass matrix M v with a 2 <-> 3 symmetry [3] can 
lead to a maximal z^-i^ mixing. Of course, a rigorous requirement of the 2 <-> 3 symmetry for 
the neutrino and charged lepton sectors cannot give a realistic lepton mixing matrix U because 
of SU(2)i symmetry [S]. We have to bring a breaking term into the 2 *-* 3 symmetric mass 
matrix by hand. A quark and lepton mass matrix model [6] has been proposed, in which the 
2-3 symmetry is broken only by a phase parameter: A mass matrix Mf is given by a form 

M f = P{5 f )M f P\5f), (1.1) 

where 

P(S) = di&g(e i5l ,e i52 ,e i53 ), (1.2) 



M f = «/ b f -c f , (1.3) 




where a/, bf and cj are real parameters. [We have defined the expression (1.3) such as > 
TOf 2 > — m fi > for bf > Cf > 0. See Eq.(A.4) in Appendix.] Here, the matrix Mf is exactly 
symmetric under the 2 <-> 3 flavor exchange, but the 2-3 symmetry for the mass matrix Mf is 



broken by the phase factor (5/3 — 5/2 7^ 0. A short review of the diagonalization of the Herimitian 
mass matrix (1.1) is given in Appendix. 

In this paper, we investigate CP violation in the lepton sector. Various models in which 
the 2-3 symmetry is broken only by phase parameters can be considered. The most naive way 
is to consider a model M u = P{5 U )M U P{8 V ) similar to the form (1.1), and such a model has 
been discussed in Ref.[6]. However, in the present paper, we take another idea that the CP 
violation in the lepton sector is cased by phase parameters which cannot be removed by a form 
My = P{5 v )M y P{5 v ) differently from those in the quark sector. The simplest realization of such 
a model is to assume a neutrino mass matrix form 



( a v a v 

V a u -c u b u e 1 ^-^ 



(1.4) 



where a u , b v and c v are real parameters. Here, the 2-3 symmetry is broken only by the phase 
parameter \ as 

M v2 2 - M v33 



M v22 + M u33 



itanx, (1-5) 



while CP is broken by the phase parameter (ft. We consider that the existence of these phase 
parameters is characteristic only in a Majorana mass matrix and we cannot consider such the 
phases in Dirac mass matrices (we assume that those are Herimitian). The mass matrix can 
lead to a nearly tribimaximal mixing when CP is maximally broken in the lepton sector. [For 
the charged lepton sector, we assume the matrix form (1.3).] 

It is hard to solve the mass matrix model (1.4) analytically. Therefore, in Sec. 2, we investi- 
gate the model (1.4) numerically. The model (1.4) has four parameters a = a u /b u , c = c u /b u , (ft 
and x f° r f° ur observable quantities sin 2 29 atm , tan 2 9 so i, | fT' 13 1 2 and R = Am 2 oZ /Am 2 tm . We will 
find that if we put a maximal CP violation hypothesis, we can obtain interesting predictions. 
Finally, Sec. 3 is devoted to summary and concluding remarks. 

2 Lepton mixing matrix 

The lepton mixing matrix U is given by 

U = U$U„, (2.1) 

where U e and U v are defined by 

U\M e XJ e = diag(m e , m M , m T ), (2.2) 

UlM u U* = di&g(m ul ,m u2 ,m u3 ). (2.3) 

A form of U e can readily be given by the analytic form U e = P{5 e )R e given in Eq.(A.l). In 
contrast to the case of U e , it is hard to express U v with an analytic form. Therefore, in this 
paper, we discuss the parameter dependences of the lepton mixing matrix U by numerical studies. 



Since, for U e , we fix the mixing angle 9 e by inputting the charged lepton masses at fj, = mz, 
the remaining free parameters for U are a v /b v , c v jb v , <j) and x- (Since the value of 9 e is very 
small, i.e. 8 e = 3.936°, the essential rotation in the lepton mixing matrix U comes from the 
neutrino sector U v .) As stated in the previous section, the phase parameter x breaks the 2 <-> 3 
symmetry, and <f> causes the CP violation. In contrast to the phase parameters (f> and x i n the 
lepton sector, there are two phase parameters a and (3 in the quark sector, which are defined 
by Eq.(A.8). The parameter (5 is only related to the breaking of the 2 3 symmetry and a is 
only related to the CP violation as shown in Eq.(A.9). However, note that those parameters in 
the quark sector do not affect quark mass values, while the phase parameters (f> and x affect not 
only the neutrino mixing matrix U u but also the neutrino mass values. In this paper, we will fix 
the phase parameters <j> and x by requiring a maximal CP violation in the lepton sector. The 
remaining parameters a v jb v and c v jb v can, in principle, be fixed by the observed values of 
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Am2 sol 

Am lt m 



m 



u2 



m ll 



m 



i>3 



m 



(2.4) 
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(we define as m 2 u2 > m 2 -^ and tan 2 # so / = |?7i2| 2 /|£/ii | 2 . Then, we can predict values of 
the rephasing invariant [7j J = Im(C/nC/22^ r i2^ r 2i)' an d the neutrino oscillation parameters 
sin 2 29 atm = 4|[/23| 2 |^33| 2 , and |C/i3| 2 . However, the observed values of R and tan 2 6 so i are not 
so accurate, at present, to determine the parameter values a = a v /b u and c = c u /b u . Therefore, 
in this paper, we will show the predictions of \ J\, sin 2 26 a t m , tan 2 6 so i, \U\^ and R as a function 
of the parameters a and c under the requirement of the maximal CP violation. 

We consider that the phenomenological mass matrix form (1.4) is given at an energy scale 
= Mz {Mz is the Z-boson mass). Since the present model gives the neutrino mass hierarchy 
m 2 x <C m 2 2 <C m 2 3 (not m 2 ^ ~ m 2 2 <C nt 2 3 ), we can regard that the renormalization group 
equation effects are not so large for the mass ratios and mixing in the present model (see, for 
example, Ref.[8j). Therefore, the numerical calculation will be done at [i = Mz, though we 
suppose that a breaking of the 2-3 symmetry is caused at a TeV energy scale. 

Note that in the following subsections A-C, we sometimes show figures with specific input 
parameter values. However, those input values are temporary ones for convenience to demon- 
strate parameter dependence of the model, and the values do not means final best-fit values in 
the present analysis. The final values are obtained from iteration of the analysis A — > B — > C. 

A. Behavior of J versus phase parameters 

Our interest is in the maximal CP violation. The maximal CP violation is defined as follow: 
the rephasing invariant J takes its maximal value for the CP violating mass matrix parameters. 
The CP violating mass matrix parameter is a phase parameter a defined by Eq.(A.8) in the 
quark sector [the model (1-1)], while it corresponds to the phase parameter <p in the present 
model (1.4). We illustrate (^-dependence of J for typical values of x i n Figs.l. Since <p ~ > ~4> 
(and x ~^ ~x) means U u — > U*, so that it means U — > U* because U e is real. Note that 
cases with (a, c) and (—a, — c) give identical results for the observable quantities, while those 
give different results for cases with (a, — c) and (— a, c). As far as the shape of |J| versus cj) is 
concerned, the cases with ac < and ac > are similar, but the maximal values of \J\ are 



different from each other as seen in Figs.l (a) and (b). Hereafter, we will refer the cases with 
ac < and ac > to as Cases A and B, respectively. As seen in Fig.l, the maximal |J| takes 
place at ~ ±13° for the cases x = 90° and \ = 120° (and 4> ~ ±12° for the case x = 60°) 
for both cases A and B. (Thus, the maximum point is slightly dependent on the parameter x-) 
The case <fi ~ 13° (0 ~ 12°) gives a maximal J with J > 0, while 4> — —13° (<fi ~ —12°) gives a 
maximal \J\ with J < 0. In this section, we take a standpoint that the maximal CP violation 
means a maximal J (J > 0). Hereafter, we confine ourselves to investigating a parameter region 
< (f) < 7r/2. Especially, since our interest is in the maximal CP violation, we will investigate 
parameter dependences for the case x = 7r /2 in which the 2-3 symmetry is maximally broken. 
On the other hand, in the limit of (j> — > 0, a case x ~~ > 7r /2 (and also x — ► — vr/2) corresponds 
to a maximal 2-3 symmetry violation as seen in Eq.(1.5). Therefore, we confine ourselves to 
investigating a parameter region < x < tt/2. 

B. (a, c) dependence 

Next, we investigate (a, c)-dependence of the observable values for sin 2 2# atm , tan 2 6>,, ;, 
1 1 2 , and R in order to give more realistic predictions of them. Since our concern is to predict 
the observable quantities at a point in which CP is maximally violated, we illustrate the (a, c) 
dependence at % = 7r /2 and 4 1 = 13°. (The value of = 13° is not a final result, because the 
value of 4> which gives a maximal \ J\ in Fig.l is dependent on the values of a and c.) 

Correspondingly to Cases A and B, a-dependences of the predicted values are illustrated 
in Figs. 2 (a) and (b), respectively, where the parameter c is fixed at a typical value c = 0.95 in 
both figures (a) and (b). We find that the present observed data [H[2] 

tan 2 9 sol = 0.469±{j;{j£ (9 aol = 34.4±i;§ degrees) (2.5) 

require a = — (0.26^q'q 2 ) and a = +(0.29^q'q3) for Cases A and B, respectively, and 

_ (7.59 ± 0.21) x 10-5 eV 2 _ . 43 2 

^ " (2.19l8;l|) x 10-3 eV 2 " (3 - 47 -°- 26) X 10 (2 ' 6) 

require a = -(0.26 ± 0.03) and a = +(0.29+^3) for Cases A and B, respectively. From overall 
view, in Case A a value a = —0.255 can give favorable predictions of tan 2 Q a tm an d R, so that 
we will use the the value a = —0.255 for Case A hereafter. For Case B, we will use a value 
a = +0.270 hereafter, although the value does not so excellent agreement with the observed 
values (2.5) and (2.6) as compared with Case A. 

Also, c-dependence is illustrated in Fig. 3 (a) and (b), which correspond to Cases A and B, 
respectively. Here, the parameter a is fixed at a = —0.255 for Case A (and a = +0.270 for Case 
B) from the results in Fig. 2. We find that the predicted values of tan 2 6 so i and R are not so 
sensitive to the value of c. If we use a parameter value c = 0.95 for both Cases A and B, which is 
favorable to the data, we obtain predictions (tan 2 9 sol , R) = (0.468,0.0347) and (0.443,0.0350) 
correspondingly to the input values (a = —0.255, c = 0.95) and (a = +0.270, c = 0.95). 

C. Phase parameter dependence 



Again, in Fig. 4, we illustrate behaviors of the predicted values J versus the phase parameter 
<fi, together with the behaviors of sin 2 29 a t m , t&n 2 6 so i, \Ui 3 \ 2 and R under the maximal 2-3 
symmetry violation x = tt/2. Here, the parameter values (a, c) = (—0.255,0.95) and (a,c) = 
(+0.270,0.95) are taken in the figures (a) and (b), respectively, from the study in the previous 
subsection. 

From Fig. 4, we find that the values of sin 2 2Q a t m and | C7" 13 1 2 are insensitive to the value of 
<p, while the predicted values of tan 2 9 so i and R are highly sensitive to the value of (f). Exactly 
speaking, in Fig. 4, since the parameters a and c have been taken so that the values tan 2 6 so i 
and R are reasonably fitted, the values of tan 2 9 so i and R are not predictions in Fig. 4. Our 
predictions are for sin 26 atm and \U 13 \ 2 : sin 2 2d atm ~ 0.98 and \U 13 \ 2 ~ 0.007 for Case A, and 
sin 2 20 atm ^ 0.93 and \U 13 \ 2 ~ 0.035 for Case B. 

D. Predictions under the maximal CP violation hypothesis 

It is interesting to assume a hypothesis of the maximal CP violation. As seen in Fig. 4, 
although the value of Jmax-, J m ax — 0.039, in Case A is larger than J m ax — 0.018 in Case B, it 
does not mean that Case A is ruled out under this hypothesis. The maximal J is required only 
for the phase parameters, not for the mass matrix parameters a v , b v and c v . Therefore, both 
Cases A and B are allowed under the maximal CP violation hypothesis. If we require a larger 
value of J, we have to choose Case B. However, the Case B predicts sin 2 29 atm < 0.93, which 
contradicts the present atmospheric neutrino data [9| sin 2 26 a t m > 0.96. Therefore, we conclude 
that Case B should be ruled out, so that the mass matrix parameters a u and c v must have the 
opposite sign to each other. In contrast to this conclusion, note that a sign of the parameter aj 
in a model (1.1) can be taken freely as seen in Eqs. (A.l) - (A. 4). 

In Fig. 5, we illustrate the more detailed behavior of the physical observable quantities under 
the maximal CP violation hypothesis. Here, although Figs. 5 (a) and (b) are similar to Fig. 2 (a) 
and Fig. 3 (b), respectively, Figs. 5 (a) and (b) are illustrated under the maximal CP violation 
hypothesis, so that the value <j) is always taken as it gives the maximal J for each parameter set of 
a and c. From Fig. 5, we conclude that sin 2 29 a t m — 0.98 and |£^i3| 2 — 0.007 for the present data 
(2.5) and (2.6). In Fig. 6, we also illustrate the contour plots for sin 2 29 a t m (Fig. (a)) , | tT" 13 1 2 x 10 
(Fig.(b)), and J (Fig.(c)) in addition to R X 10 and tan 2 8 so i given in Eqs. (2.5) and (2.6) in the 
(a, c) parameter plane under the maximal CP violation hypothesis for each parameter set of a 
and c. 

4 Summary 

In conclusion, we have assumed a neutrino mass matrix in which a 2-3 flavor symmetry 
is violated only by a phase parameter, and which can lead to a nearly tribimaximal neutrino 
mixing, and thereby, we have investigated behaviors of the neutrino oscillation parameters under 
the maximal CP violation hypothesis. We find that the predicted values of sin 2 29 a t m and | 3 1 2 
are insensitive to the input values R and tan 2 6 so i for R ~ 0.035 and tan 2 6 so i ~ 0.5: 

sin 2 29 atm ~ 0.98, \U 13 \ 2 ~ 0.01. (4.1) 

Especially, we are interested in the case with \R\ = 0.0321 which is a lower bound of the observed 
value of \R\ as shown in Eq.(2.6), because the value gives us a upper bound of the predicted 



value for sin 2 26 atm and a lower bound of the predicted value for \Uis | 2 . As seen in Fig.6, our 
model rules out the following values 

sin 2 2# atm > 0.989, | 1 2 < 0.0046, (4.2) 

if neutrino oscillation data establish \R\ > 0.321. On the other hand, if the present upper bound 
of \R\, \R\ = 0.0390, is established, the following regions 

sin 2 29 atm < 0.978, \U 13 \ 2 > 0.0098, (4.3) 

are also ruled out. Those bounds (4.2) and (4.3) are within reach of near future experiments. 
Therefore, if future experiments give us values in the regions (4.2) and (4.3), the present model 
will be ruled out. On the contrary, if future experiments report observed values within 

0.989 > sin 2 26 atm > 0.978, 0.0046 < \U 13 \ 2 < 0.0098, 

0.021 > J > 0.013, (4.4) 

correspondingly to 0.0321 < \R\ < 0.0390 and 0.428 < tan 2 6 sol < 0.516, the neutrino mass 
matrix form (1.4) will become a promising candidate of the neutrino mass matrix form. 

Also, note that the predicted value of tan 2 9 so i is highly correlated to the predicted value of 
R under the values (4.1). This result is shown in Fig.7 by using only the observable quantities 
with fixing typical values for sin 2 26 atm such as sin 2 29 atm = 0.980 and 0.985. This figure is also 
useful as a touchstone of the present model. 

The neutrino mass matrix (1.4) is a very simple form and it has phase parameters which 
are characteristic only in a Majorana mass matrix. We expect that the results will be checked 
in the very near future experiments. 
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Appendix 

We give a short review of the diagonalization of the mass matrix (1.1) and the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix as an example of the flavor mixing matrix. 
(Although we have denoted as (M/)n = for simplicity, it is not essential for the following 
discussions. For the case (M/)n / in Eq.(l.l), it is only necessary that the parameter values 
bf and 6/3—6/2 are redefined under a common shift of the mass eigenvalues m/i — > mfi+(Mf)u.) 

The matrix Mf given in Eq.(1.3) is diagonalized by a rotation 



R f = R 1 (~)R 3 (d f ) 



( cos 9/ s'm9f \ 
"72 sin ^ 75 COS0 / "75 



(Al) 



as 

RjMfRj = Df = diag(m/i,m/ 2 ,m/ 3 ), 



(A.2) 



where 



COS Uf 



m f2 
m/ 2 — m/i 



sin ( 



/ —mfi 
m/2 - 771/1 



(A3) 



m n = hi b f ~ c f ~ \/ 8a f + ( b f ~ c /) 2 » 



m f2 = \{hf-c f + ^8a 2 f + (6/ - c/) 2 J , 
m f3 = bf + c f . 



(A4) 



Here and hereafter, we use the following notations for the rotation matrices 
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, #2(0) = 
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Vo 


—scj 






c / 




V o 


1 / 



(A5) 



where c = cos# and s = sin#. When we apply this model to quark mass matrices, we obtain 
the CKM quark mixing matrix V as follows: 



V = UlU d = R T u P{5)R d = Rj(e u )Rj(~Tr)P(S)R 1 (~7r)R 3 (e d ), 



(A6) 



where 5i = S u i — 5di (i = 1,2,3) and tan#/ = y /— m fi/ m f2- We can rewrite the factor 
rT(-\v)P{5)Ri{-\-k) as 



1 



1 



R[(— 7T )P(5)R 1 {--tt) 



( e iSl 

\ (e iS2 + e iS3 ) \ (-e iS2 + e i<53 ) 
\ \ (-e i52 + e iSs ) \ (e i52 + e i5s ) 



( e~ ia 

cos/3 isin/3 | = e l{a+5 ^P l (-a)P 3 (- l -^)R l (f3)P 3 {]^), (A7) 



\ i sin (5 cos /3 



where 



a = ^(S 3 + S 2 ) - 5 U (3 



(A 



Pi (5) = diag(e i<5 , 1, 1) and P 3 (6) = diag(l, 1, e iS ). Since P 3 and R 3 are commutable each other, 
the phase matrices P 3 (— 7r/2) and P3(7r/2) can be eliminated by redefining phases of the up- and 
down-quarks, respectively, so that we can obtain a CKM matrix form as follows: 



V = R 3 r (e u )P 1 (-a)R 1 (P)R 



3\yd) 



(A9) 



This expression (A. 9) is well-known as the Fritzsch-Xing expression [10J of the CKM matrix. 
Note that the phase parameter a still plays a role as a CP violating phase, while the parameter 



plays only a role as a rotation parameter 623. It is known that the model is in favor of the 
observed values of the CKM parameters under a hypothesis that the CP violation is maximally 
violated, i.e. a = n/2. 
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Fig. 1 Behavior of the rephasing invariant J versus the CP violation parameter 
4>. Figures (a) and (b) are illustrated for typical values (a, c) = (—0.27, +0.95) 
and (a, c) = (+0.27, +0.95), respectively. The dashed, solid, and dotted curves 
correspond to cases with \ = 120°, \ = 90°) an d X = 60°, respectively. 




-0.35 -0.33 -0.31 -0.29 -0.27 -0.25 -0.23 -0.21 

a 



1.1 
1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



(a) 






sin 2 2 a tm 




tan Osol 


- io|u 13 | 2 

1 I 1 1 


10R 



0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34 



a 



(b) 



Fig. 2 a-dependence of the predicted values for sin 29 atm , tan 6 so i, R x 10 and 
j 1 2 x 10. Other parameters are fixed at typical values c = c v jb v = 0.95, x = 7r /2, 
and (f) = 13°. Figures (a) and (b) correspond to Case A with ac < and Case B 
with ac > 0, respectively. 
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Fig. 3 c-depencence of the predicted values for sin 2 2# aim , tan 2 6 so i, R x 10, and 
\U 13 \ 2 x 10. Fi gures (a) and (b) correspond to Case A with ac > and Case B with 
ac < 0, respectively and those are illustrated by taking input values a = —0.255 and 
a = +0.27, respectively. Other parameters are fixed at typical values % = ^/^ an d 
6= 13°. 
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Fig. 4 Behavior of J, sin 2 20 atm , tan 2 9 so [, R x 10, and \U\^\ 2 x 10, versus the CP 
violating phase parameter <j) under the maximal 2-3 symmetry violation \ = it/2. 
Figures (a) and (b) correspond to Case A with ac < and Case B with ac > 0, 
respectively, and those are illustrated by taking input values (a,c) = (—0.255,0.95) 
and (a, c) = (+0.27,0.95), respectively. 
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Fig. 5 o-dependence (Fig. (a)) and c-dependence (Fig. (b)) of the predicted values 
for J, sin 2 29 a tm, tan 2 9 so i, R x 10, and |£7i3| 2 X 10 under the maximal CP violating 
hypothesis: the value of <f> is taken as J is maximal for each parameter set of a and 
c. Here, the phase parameter x is fixed at x = 90° a t which J takes the maximal 
value, and the parameter c is taken as c = 0.95 for Fig. (a) and a as a = —0.255 for 
Fig. (b). 
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Fig. 6 Contour plots of predicted values of sin 2 26 a t m (Fig. (a)) , \Uis\ 2 x 10 
(Fig. (b)), and J (Fig. (c)) in the a-c parameter plane under the maximal CP 
violating hypothesis: the value of (j> is taken as J is maximal for each parameter set of 
a and c. Here, the phase parameter \ is fixed at % = 90° a t which J takes the maximal 
value. We also present contour plots of R x 10 and tan 2 8 so i which are correspond to 
the experimental constraints, 0.0321 < \R\ < 0.0390 and 0.428 < tan 2 6 sol < 0.516, 
respectively. Here the darker region corresponds to smaller contour values. Note 
that the overlapped region of the contours of WR and tan 2 9 so i is the allowed region 
of the a-c parameter plane of the model. 




Fig. 7 Behavior of tan 2 9 so i vs R under the maximal CP violation hypothesis 
with taking typical values of sin 2 29 a t m . The solid and dashed curves correspond to 
the cases with sin 2 29 at m = 0.980 and sin 2 29 atm = 0.985, respectively. The predicted 
values of \U 13 \ 2 for sin 2 29 atm = 0.980 and sin 2 29 atm = 0.985 are \U 13 \ 2 ~ 0.0087 and 
|^13| 2 — 0.0066, respectively, in this parameter region of (a,c). 



